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We study Nf = 2 lattice QCD with nonperturbatively improved Wilson fermions at finite temperature on 16^ • 8 
lattices. We determine the transition temperature at ^ ~ 0.8 and lattice spacing as small as a ~ 0.12fm. The 
string breaking at T < Tc is also studied. We find that the static potential can be fitted by a simple expression 
involving string model potential at finite temperature. 



1. Introduction 

Recent studies of TV/ = 2 lattice QCD at fi- 
nite temperature with improved actions have pro- 
vided consistent estimates of Tc- Bielefeld group 
employed improved staggered fermions and im- 
proved gauge field action CP-PACS collabo- 
ration used improved Wilson fermions with mean 
field improved Csw and improved gauge field ac- 
tion 1^. Both groups were able to estimate Tc 
in the chiral limit and their values are in a good 
agreement. Still there are many sources of sys- 
tematic uncertainties and new computations of Tc 
with different actions are useful as an additional 
check. We made first large scale simulations 
of the nonperturbatively 0(a) improved Wilson 
fermion action. Moreover we performed simula- 
tions with the lattice spacing a much smaller than 
in studies by Bielefeld and CP-PACS groups. Our 
small lattice spacing helps us to determine param- 
eters of the static potential in fidl lattice QCD at 
finite temperature. Our other goal is to study 
the vacuum structure of the full QCD at T > 0, 
in particular a relation between string breaking 
observed with the help of the Polyakov loop cor- 



relator at T < Tc and abelian monopoles. 

The fermionic action we employ is of the form 



Sp — S 



(0) 



where S*^' is the original Wilson action, Csw is 
determined nonperturbatively ||^ . We use Wilson 
gauge field action. 

We took advantage of availability of T = re- 
sults obtained with the action (|l]) by UKQCD 
and QCDSF collaborations |^ to fix the physical 
scale and — ratio. Studies made by these coUab- 
orations also confirmed that 0(a) lattice artifacts 
are suppressed as expected [|| . We then may hope 
that lattice discretization errors of our results are 
small. So far only iVj = 4 and 6 finite temperature 
results obtained with this action are available . 
These results were obtained at rather large quark 
mass (rnT^/nip > 0.85) and lattice spacing. 

Since Csw is only known nonperturbatively for 
/3 > 5.2, we have to work at rather large values 
of (3. Thus to study transition at smaller values 
of the quark mass we need to go to larger Nt 
values. This dictates our choice Nt = 8. We 
choose spatial extension of the lattice Ns = 16 
as a compromise between computational burden 
and need to reduce finite size effects. 



*Talks given by V. Bornyakov and Y. Nakamura at Lattice 
2002, MIT, Cambridge MA, USA. 
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Table 1. Simulation statistics. 

It is known that in quenched QCD the order 
parameter of the finite temperature phase transi- 
tion is Polyakov loop and corresponding symme- 
try is global Z{3) symmetry. In chiral QCD the 
order parameter of the chiral symmetry breaking 
transition is the chiral condensate (ipip). There is 
no phase transition at intermediate values of the 
quark mass, only crossover. As numerical results 
show 0] both order parameters can be used to lo- 
cate the transition point. We use only Polyakov 
loop leaving computation of the chiral condensate 
for future studies. 

2. Simulation details 

We use Hybrid Monte Carlo algorithm with pa- 
rameters 5t — 0.0125, rir = 20, with acceptance 
rate about 70%. Simulations were done on Hi- 
tachi SR8000 at KEK and MVS lOOOM at RAS, 
Moscow. We needed from 1000 (r = 250) to 3000 
(r = 750) trajectories for thcrmalization, depend- 
ing on K and Our simulations were performed 
for /? = 5.2, 5.25. The values of k and correspond- 
ing number of trajectories are listed in Table 1. 

We measured the Polyakov loop and plaque- 
tte on every trajectory. They were used to mea- 
sure corresponding susceptibilities. Depending on 
K every 5th, every 10th or every 20th trajectory 
were used to compute Polyakov loops correlator 
{LsLI), where Ls = TiPs, Ps = Ila^o ^0;xo,x- 
The MA gauge was also fixed on these configu- 
rations. Simulated annealing algorithm has been 
used to fix the gauge M. Abelian and monopole 



A 
V 




0.8 r 
0.7 - 



A 
V 



o _^ mou 

' ~ 1 ^ f^p'" 

. 5 3^mon X a^ph 



>-+--: 



0.4 
0.3 - 



+ HHH- + 



Figure 1. Nonabelian Polyakov loop at both /?'s 
(top), and the abelian, monopole and photon 
Polyakov loops at /3 = 5.2 (bottom). 

Polyakov loops, their susceptibilities and correla- 
tors were measured on gauge fixed configurations. 
It turns out that for abelian and monopole ob- 
servables the signal/noise ratio was better than 
that for gauge invariant nonabelian observables. 
This observation is in agreement with results from 
quenched QCD at T > and both quenched and 
unquenched QCD at T = 0. For this reason we 
use abelian and monopole observables to deter- 
mine transition temperature and to study string 
breaking. To obtain more precise results with 
our limited statistics we also employed hypercu- 
bic blocking introduced recently 

3. Transition temperature 

In Fig.|l| we show results for average of various 
kinds of Polyakov loops. One can see that (L) 
is a smooth function of k. In the bottom part 
of this figure we show that as in quenched the- 
ory (LAb) ~ (imon)(iph) and that (Lph) does 
not show any changes at the transition. We de- 
termined Kt, the critical k, from the maximum 
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Figure 2. Nonabelian Polyakov loop susceptibil- 
ities (top), and abelian and monopole Polyakov 
loop susceptibilities at /3 = 5.2 (bottom). 

of the Polyakov loop susceptibility. We found 
Kt = 0.1344(1) at /3 = 5.2 and K( = 0.1340(1) 
at /3 = 5.25, see Fig. § 

An interpolation formula for ro/a Q gives us 
Tcro = 0.54(2) and 0.56(2). Or, using r^^ = 394 
MeV, we obtain the critical temperature in phys- 
ical units, = 213(10) and 220(5)MeV. Us- 
ing again the data from |^ we get mTr/wp = 
0.78,0.82 at the transition points. The suscepti- 
bilities for Abelian and monopole Polyakov loops 
have maxima at the same k, see Fig. 2. 

In Fig. ^ our results for transition temperature 
are shown in comparison with those of Refs. [Q, 
|. Thus our results are in a good qualitative 
agreement with the results of the Bielefeld group. 

4. String breaking 
4.1. T = 

At r = string breaking has been observed 
so far with mixed operators including explicitly 
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Figure 3. Transition temperature. 

static-light meson state and no sign of string 
breaking from Wilson loop has been found [|lOj. 
The last fact is believed to be due to very poor 
overlap of the broken string state with the Wilson 
loop operator. Keeping two terms in the Wilson 
loop spectral representation, corresponding to the 
string and two-meson states, one gets 

Wir, t) = Cy(r)e-(^'+^="-M)* + CE{r)e~''^-\ (2) 

where Vstrir) the usual confining potential and 
E is the static-light meson energy. The overlap 
Cy(r) is of 0(1) while Csir) might be small. A 
quantitative estimate of this overlap has been sug- 
gested in ||ll| and then derived in abelian projec- 
tion approach [Q: Csir) e"^^ ''. We believe 
that this is correct estimate at least in strong 
coupling, small n limit. If we assume that the 
real overlap is indeed of this order then the ef- 
fect of the second term in eq. (^ would become 
detectable for both r and t of order of 2 fm, 
i.e. for very large Wilson loops with very small 
value. The real string breaking distance Vsb is 
to be estimated from the equality of two expo- 
nents in eq. (||): 2m — a ■ rgb — Tr/{12rsb), where 
m = E — Vo/2 may be called the string breaking 
energy or the effective quark mass ||l^ . The Wup- 
pertal group result is Vsb = 2.3ro at m.„/mp ~ 0.7 
im and the CP-PACS result is r^b = 2.2ro at 
rriTr/mp = 0.6 ^t| . Taking the string tension in 
full QCD from Ref. 0, y/^ro = 1.15(1), we get 
the effective quark mass 2m ~ 2.9/ro ~ l.lGeV. 
This is in a surprisingly good agreement with es- 
timate made in |T3] using a different approach. 
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4.2. T > 

At non-zero temperatures the string breaking 
has been observed in The heavy quark po- 
tential V is related to the Polyakov loop cor- 



relator, e" 



-V{r,T)/T 



{LsLi)/9. 



In the limit 



|a; — yl oo, (LgLt) approaches the cluster value 
|(L)P, where |(L)p ^ because the global Z3 
symmetry is broken by the fermions. 
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Figure 4. Heavy quark potential from nonabelian 
(top) and monopole (bottom) Polyakov loops, fit- 
ted by eq.(||) and eq.(|^), respectively. Horizontal 
lines show — 2roTlog(L/3) with the error bars. 

The spectral representation for the Polyakov 



loop correlator is |16 



much higher energies one can write 

1 
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-V{r,T)/T 



Va.t{r,T)/T 



In full QCD wc must take into account the ex- 
cited states. As it was discussed above in full 
QCD at T = the singlet potential can be de- 
scribed by string model potential up to some dis- 
tance riir- Beyond this distance another state 
becomes the ground state of the system of two 
static quarks which one can call broken string 
state or two heavy-light mesons state. So there 
are two distinct states in the spectrum, one of 
which becomes the ground state at proper dis- 
tance. The situation must be similar at small 
temperatures. We now assume that at temper- 
atures T < Tc the Polyakov loop correlator can 
be described with the help of these two states, 
namely string state and broken string (two me- 
son) state. In this paper we do not single out the 
singlet potential. We understand that the contri- 
bution from the octet potential may contaminate 
our results. The work on calculation of the sin- 
glet and octet potential separately is in progress. 
We now adopt the following representation for the 
Polyakov loop correlator 



1 



'{Vo + V,tr{r,T))/T 



-2E{T)/T 



(3) 
(4) 



g \ 1,/ 

Vstr{r,T) = -(^^- iarctan(2rT) 

o-(r) H arctan r H In 1 

^ ^ 3 2rT) 2 V 

E{T) = Vb/2 -f m{T) , (5) 

where m{T) is effective quark mass at finite tem- 
perature. The T ^ string potential (^ was cal- 
culated in Ref . jlTj . An alternative way to fit the 
static potential in full QCD at finite temperature 
was proposed long ago [p8|: 



n=0 



with integer Wn- At T = we get V{r,T) = 
£'o(r). In contrast, V{r,T) at T > gets con- 
tributions from all possible states. In pure gauge 
theory these are singlet and octet states and their 
excitations. Assuming that the excitations have 



V{r,T) = -(1 
A* 



--) - %~ 
r 



fir 



(6) 



with parameters a and /i. In (^ only /i or 
only (7 |19| are temperature dependent. We used 
function (|6|) to fit our data and compare with our 
model eq. (|)-(|). 

Fig. U shows the nonabelian static potential for 
some values of T/Tc- The fit curves are two-state 
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Figure 5. Best fit parameters for nonabelian and monopole potentials as functions of temperature. Solid 
lines on left figure show quenched results M . 



fits introduced above. We find good agreement of 
the fit with lattice data. As it has been already 
mentioned above to reduce statistical errors we 
applied hypercubic blocking (HCB) We used 
the same parameters of HCB procedure as in [|j . 
HCB helped to reduce the statistical errors by 
about factor 2 and to improve the rotational in- 
variance. The potentials without and after HCB 
coincide within error bars up to a constant at all 
distances r apart from r — a and 

Parameters of our fit (^) are presented in Fig. |[ 
Our values for the ratio ct(T)/ct(0) are higher 
than those obtained in quenched QCD in Ref . . 
Our values for the quark effective mass are also 
higher than those obtained in Having pa- 

rameters of the potential determined we can cal- 
culate the string breaking distance r^b using the 
relation Vstrirsb,T) = 2to(T). From Fig. || one 
can see that r^f, decreases down to values ^ 0.3 
fm when temperature approaches critical value. 
Our fit using Vstr is probably not valid when 
rsb becomes so small. It still can be valid for 
T/Tc < 0.95 when r^t > 0.5/to. 

The quality of our data for (LgLt) does not 



allow us to distinguish between the fits eq. (|§|5|) 
and (|^). However, the effective string tension in 
the fit (H) is unphysically large, otq ~ 6. 

We calculated also the potential after abelian 
projection. It is known from quenched QCD stud- 
ies that the monopole part of the abelian gauge 
field gives a potential which is linear down to very 
small distances, i.e. the Coulomb coefficient is 
compatible with zero. We observed similar be- 



haviour at small distances as is seen from Fig. |j. 
At large distances the monopole potential shows 
screening as the nonabelian potential does. We 
fit the potential similarly to eq.(||): 

1 mon 

-{LsLi) = e-(^o"" + K"r('-.T))/T^g-2£— (T)/T(7) 

9 ^ 

= o-'^^r and i;'°°"(T) = VJ"™ -l-m'°°"(T). 
The fit is shown in Fig.^ by the solid line. 



^^ • Full QCD 

- Quenched QCD 



Figure 6. Density of magnetic currents. 



5. Monopole Dynamics 

We studied monopole density p and asymmetry 
in the monopole density i] = {pt — ps)/{pt + Ps) 
where pt{ps) is density of time-like (space-like) 
magnetic currents. In Fig. ^ we compare the 
monopole density in full and quenched QCD. Our 
results for quenched QCD were obtained on 16^-8 
lattice. The density in full QCD is substantially 
higher than that in quenched theory in agreement 
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with our earlier results at T = |2C[|. The be- 
haviour of T] is qualitatively similar in both the- 
ories, i.e. r] is close to zero below transition and 
increases with temperature above Tc- 

6. Conclusions and Outlook 

Our results for Tc obtained on 16^ • 8 lattice 
at ^ ~ 0.8 and lattice spacing small in compar- 
ison with work of Bielefeld group are in agreement 
with their results |^ . To test finite size effects and 
to determine transition temperature closer to the 
chiral limit we are planning to make simulations 
on • 8 and 24^ • 10 lattices. 

Heavy quark potential has been measured in 
both phases. We introduced two-state param- 
eterization eq.(|3[^) for the heavy quark poten- 
tial in confinement phase and found good agree- 
ment with numerical data obtained in that phase. 
Using this parameterization we computed string 
tension, quark effective mass and string breaking 
distance. We found that the ratio a{T)/a{0) de- 
creases toward the value 0.7 when T approaches 
Tc- This value is substantially higher than the 
value obtained in quenched QCD. Whether this 
is due to systematic effects (meson-meson inter- 
action, contribution of the octet potential, and 
some others) or is a physical effect deserves fur- 
ther study. Our results for the quark effective 
mass show good qualitative agreement with ear- 
lier results [Q. The results for string breaking 
distance imply that our fit might be inappropri- 
ate for temperatures T/Tc > 0.95. 

We observed string breaking at T < Tc also 
from abelian and monopole Polyakov loop corre- 
lators. We found that string tension, quark ef- 
fective mass and string breaking distance calcu- 
lated from the monopole Polyakov loops correla- 
tors are in agreement with the values obtained 
from the nonabelian correlators. Thus we found 
abelian and monopole dominance in the confine- 
ment phase of full QCD. 
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